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A study of colour field distributions in the baryon 
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The distributions of chromo-electric and chromo-magnetic field associated with flux tubes in the baryon are 
studied in SU(3) lattice QCD. Maximal Abelian projection is used to reduce the statistical fluctuations. For 
a fixed source geometry, many different string configurations are possible. We investigated whether the string 
configuration, that is the choice of operator, biases the observed flux distribution. 



1. Introduction 

Recently, the problem of how chromo-electric 
and chromo-magnetic fields are distributed in a 
baryon has been of considerable interest. It has 
been investigated in terms of the potential |1I2| . 
But the flux distribution can be calculated di- 
rectly using lattice QCD. There are many studies 
for a meson [31415) . but only a few for a baryon. 
The first attempt to do this was by Flower[B] but 
it was very noisy and not conclusive. However, in 
the past year, Ichie et al.]J\ have presented very 
nice results using Abelian projected fields. 

Generally speaking, all lattice simulations have 
some systematic errors and it is important to have 
some estimate of these errors in order to assess the 
results. What are the sources of systematic error 
in the map of flux distribution? 

One is interested in the shape of the flux distri- 
bution. However, in order to do the calculation 
some choice for the three quark operator has to 
be made. This operator consists of three quarks 
connected by strings of gauge field links which as- 
sume some shape. The systematic effect which we 
study here is whether the shape of the operator 
influences the measured flux distribution. 

Of course, the ground state property in a lat- 
tice simulation should be independent of the op- 
erator. But the conditions needed to obtain this 
ideal situation may be different for different quan- 
tities and may be difficult to achieve in practice. 



This can be seen by considering the time depen- 
dence of the two-point and three-point correlation 
functions. 

For the two-point function 
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For the three-point function (T > t) 
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The generic form of the correlation function 
which gives the field distribution is 
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In this case, the condition T >> t >> is re- 
quired to isolate the ground state contribution. 
For the two point function, which yields the po- 
tential, one only needs T >> 0. If a simulation 
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Figure 1. The three quark Wilson loop (a) T- 
shape string configuration (b)L-shape string con- 
figuration. 



The simulation is done in quenched QCD using 
the Wilson plaquette action with /3 = 6.2 on 24 4 
lattices. In order to be able to study systematic 
effects it is necessary to control statistical errors. 
The multi-hit procedure was used for temporal 
links in order to reduce the statistical fluctuations 
in the potential calculation. In addition, smear- 
ing of spatial links is used to enhance the ground 
state. The flux distribution calculation is very 
noisy. Ichie et al.\J\ have shown that maximally 
Abelian projection, which preserves the correct 
long distance behavior of the Wilson loop, can be 
used very effectively to investigate the flux tube 
in the baryon. We adopt this method here. 

Fig. 1 shows the choice of the operator. Three 
quarks are located on the x — y plane with the 
distance k from the origin, {Ql(x, y), Q2, Q3) — 

((-M),(0,*),(M)) m (a), ((0,0),(0,fc),(fc,0)) 
in (b). The operator (b) has not been considered 
before. 

3. Results and Discussion 



shows a dependence on the choice of operator this 
indicates that there are non-ground state contri- 
butions. 

2. Simulation 

We studied the three quark potential and the 
field distributions. The three quark Wilson loop 
is defined as 

W 3Q = y abc t alblc ,ufufuf (4) 
with the path-ordered link variables 
Uj^Pexpfcgf dx^(x)} (j = 1,2,3). (5) 

The path is denoted by Tj in Fig. 1. 

The correlation function can be described as 



(P^(x,t)W 3Q (T)) 



(P^(t)). (6) 



Here, P M „ denotes an unsmeared plaquette, used 
as a probe of the chromo-electric or chromo- 
magnetic field. In our simulation the plaquette 
insertions were made at t = T/2. 
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Figure 2. The three quark potential. 



Fig. 2 shows the three quark potential for dif- 
ferent operators. L min describes the total min- 
imum length from the physical junction to each 
quark. The circles denote the case of three quarks 
located on the axes at the same distance k from 
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the origin, while the triangles and the squares 
show the quark geometries of Fig.l (a) and (b) 
respectively. There is no evidence for any opera- 
tor dependence. 
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Figure 3. The action density (a) T-shape string 
configuration (b) L-shape string configuration. 



Fig. 3 shows the action density for operators 
with different string configurations. The three 
quark positions, ((— 6, 0), (0, 6), (6, 0)) in (a) and 
((8,0), (0,0), (0,8)) in (b) are chosen since the 
inter-quark separations are similar. T = 8 is 
considered since it is large enough to extract the 
ground state for the potential. There a visible 
difference between the shapes of the distributions 
which are still biased by the form of the operator. 
This operator dependence suggests that excited 



states are still contributing. 

Ideally, the physical ground state property in 
a lattice simulation should be independent of the 
choice of the operator. For the three quark poten- 
tial, this seems to be readily achievable. However, 
since the flux distribution contains the three- 
point function, it requires large time separation 
not just between the quark source and sink but 
also between the source and sink and the flux 
probe. These conditions are not easy to meet 
and in our simulations, operator dependence in 
the flux distribution is visible. It is important for 
future calculations to check for this systematic ef- 
fect, especially if it is claimed that the measured 
flux distribution has the same shape as the strings 
contained in the baryonic operator. 

To reduce the systematic error discussed here, 
even more effective techniques to reduce statis- 
tical fluctuations would be helpful. This would 
allow calculations at larger time separations. Al- 
ternatively, operators without strings, as used, for 
example, by the MILC Collaboration 8 to calcu- 
late the quark-antiquark potential, might be use- 
ful. 
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